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ABSTRACT 


In this paper a series of fracture problems in composite 
materials are identified, their methods of s°lutn° n ar e br-i ef l y 
discussed, and some sample results are presented The main 
problem of interest is the determination of the stress state in 
the neighborhood of localized imperfections such as cracks and 
inclusions which may exist in the composite. Particul ar emphasi s 
is placed on the evaluation of quantities such as the stress 
intensity factors, the power of the stress singularity, and th 
strain energy release rate, which may be used directly or indi- 
rectly in connection with an appropriate fracture criterion for 
the prediction of fracture initiation and propagation load _ 
levels. The topics discussed in the paper include a crack in 
layered composites, a crack terminating at and going through a 
bi-material interface, a penny-shaped crack in a filament- 
reinforced elastic matrix, and inclusion problems in bonded 

materials. 


The research reported in this work was supported by The National 
Science Foundation under Grant GK-11977 and by the National 
Aeronautics and Space Administration under Grant NGR-39-007-01 1 . 



1 . INTRODUCTION 


Generally the fracture of composite materials may be con- 
sidered from two different points of view. In the first approach 
the main interest is in estimating or studying the "bulk strength' 
of a given structure or of a part under a known system of external 
loads and environmental conditions. In this type of studies it 
is usually assumed that the material is statistically homogeneous 
and existing imperfections (which are unavoidable) are randomly 
distributed. Here, the very nature of the problem requires that 
some kind of statistical strength theory be used as a guide in 
the investigations. 

In the second approach to the fracture studies in composites, 
one is basically interested in the initiation of fracture propa- 
gation from the "localized" imperfections which are known (or 
assumed) to exist in the material. Here, before embarking on 
any elaborate analysis, for this local fracture initiation one 
has to adopt a proper fracture criterion and decide on the type 
of "load factor" to be evaluated. Usually the fracture criterion 
consists of a simple comparison between a calculated load factor 
and a material constant which is determined from certain standard 
experiments. Hence, the approach is, of necessity, deterministic 

In this paper we are interested only in the second type of 
problems, namely, the fracture problems arising from isolated 
imperfections. The most common forms of these imperfections 
which may be found in composites are voids, cracks, and inclu- 
sions with varying degrees of relative stiffness. In studying 
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the composites which consist of perfectly bonded multiphase 
elastic continua, it is usually assumed that the imperfection 
located near to or on an interface may be idealized as a crack 
or as a flat inclusion, the reason being partly- analytical . As 
in most analytical work on the subject, in this paper we will 
restrict our attention to the linearly elastic materials and 
concentrate our numerical effort on the evaluation of the strength 
of the stress singularities known as the stress intensity factors. 
These factors are defined in terms of (or, are obtained from) 
the asymptotic expressions for the stress state around the singu- 
lar points, which are generally of the following form: 


° 1 j< r > e > 


-ii_ f..V,e) ♦ f,, 2 (r.e> * 0 (r'-“) (1) 

(2r)° 1J (2r)° 1J 


where k-, and k 2 are the stress intensity factors, (r,0) are the 
polar coordinates in the plane perpendicular to the periphery of 
the crack or of the inclusion, f i:j k , (k=l,2) are bounded functions, 
and a is the power of the singularity (0 < a < 1). The conjecture 
here is that, regardless of the nature of the singularity at the 
imperfection front, in at least ideally brittle materials k ] and 
k 2 are a reasonable measure of the severity of deformations and 
stresses at the expected location of the fracture initiation, and 
hence should provide sufficient information for the calculation 
of the "load factor" which may be used in a proper fracture 

criterion. 


Figure 1 shows various possible modes of crack initiation 
in the neighborhood of a singular point. Figure la describes the 
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most probable crack propagation mode for a crack imbedded in a 
symmetrically loaded homogeneous (and from the viewpoint of frac- 
ture resistance, isotropic) medium. Here, as the crack propa- 
gates, the stress state around the crack tip remains autonomous 
or self-similar, the only change being in the stress intensity 
factor. Therefore, in this case, if one prefers to use it, an 
energy balance type of fracture criterion is applicable. The 
only other fracture mode for which this could be said is the 
propagation of an interface crack shown in Figure lc. Here, 

( k] 2 + k 2 2 ) is the measure of the strain energy release rate and 
may be considered as an appropriate load factor [1]. The type of 
fracture shown in Figure lb may initiate from a crack tip around 
which the stress state is not symmetric. In bonded materials 
containing an interface flaw, if the adhesive strength of the 
bond is greater than the rupture strength of one of the adjacent 
materials, the fracture mode may be that described by Figure Id. 
Figure le shows the fracture initiation around the end point of 
a flat elastic inclusion. The common feature of the stress singu- 
larities in fracture problems described by Figures la-le is that 
the power of the singularity a shown in (1) is 1/2. 

On the other hand in the crack initiation problems for the 
flaws running into and that going through a bi-material interface 
shown in Figures lf-li this no longer is the case, i.e., if 
f y 2 then a t 1/2 where y-, and y 2 are the shear moduli 

(1 > a > 1/2 for y 1 > u 2 and 1/2 > “ > 0 for U 1 < u 2 in FigureS 
lf-lg, and 1/2 > a > 0 for Figure li). In the problem of the 

flaw running into an interface, whether further fracture 
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propagation would be the cleavage of the adjacent medium ( f ) s 
debonding along the interface (g), or reflection back into the 
first medium (h) may again depend on the relative strength-to- 
Toad factor ratios for various possible fracture modes. The same 
principle also applies to the initiation of debonding crack shown 
in ( i ) . 

In this paper we will give the results regarding a, k n , and 
f-- n , ( i , j , n = 1,2) for various problems arising from the fracture 
of composite materials. As to how to use these results in prac 
tice, it may depend on one's preference for a particular fracture 
criterion. However, a certain comparative test common to all 
fracture modes shown in Figure 1 is worth noting. Whatever the 
nature of the conjectured fracture mechanism (namely, cleavage 
starting from the singular point, or fracture initiation and 
propagation through hole or craze formation and coalescence around 
the singular point), the intensity of the (tensile) stress field 
within the immediate neighborhood of the singular point will be 
the major factor controlling the fracture initiation, and this 
stress field intensity is clearly dependent on k^/r , (i - 1,2) 
(see equation (1)). Hence, comparing, for example, the maximum 
value of the cleavage stress a 0Q (in the homogeneous medium) or 
an interface stress vector computed in terms of k^ and a at a 
characteristic distance from the singular point, with the corre- 
sponding quantities obtained from idealized test configurations 
under rupture conditions may provide a rather simple criterion 
for fracture initiation. 
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2 . 


GENERAL METHOD OF SOLUTION 


/ 


Assuming that the solution to the elastostatic problem for 
the imperfection-free nonhomogeneous solid is available, the 
singular behavior of the solution around an isolated imperfection 
may be obtained from the stress disturbance or perturbation prob- 
lem. If desired, the complete solution may then be found through 
superposition. The perturbation problem is a mixed boundary 
value problem and can always be formulated in terms of a system 
of singular integral equations. In the case of plane strain, 
the generalized plane stress, and the axisymmetric problems for 
the type of imperfections shown in Figure 1 (and for many other 
problems), the general form of these integral equations may be 
expressedasfollows: 


l a-H 


i J 


*j (x) 


+ ~ 
7T 


a 

f 


-a 


f b u *j (t) 


dt 

t-x 




3 

I 

k = l 
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i jk 


*j (t) 


(x+a) k ~ 1 + 

( t+x + 2a) k 



*j (t) 
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k-1 


( t+x-2a ) k 


]dt 


a 

► 

+ 

-a 


k i j ( x , t ) 4>j(t)dt = Y i 


f^x) , 


( i =1 ,2) , -a<x<l 


( 2 ) 


where a is the half-length or the radius of the imperfection and 
a i j , b ij5 c ijk , d ijk , and Yi (U = 1.2, k = 1,2,3) are constants 
which are dependent on the elastic properties of the materials 
forming the composite. In crack problems the input functions f ^ , 
(i = 1,2) are the components of the traction vector on the crack 
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( 3 ) 


surface and the unknown functions 4> - are defined by [2-8] 

<M X > - h < u i + - u i‘> • (i ' 1>2) 

where u-, and u 2 are the crack surface displacements parallel and 
perpendicular to the plane of the crack, respectively, x is the 
(rectangular or radius) coordinate in the plane of the crack, 
and superscripts + and - refer to the opposite faces of the crack. 
Note that in this case 4> ^ (i = 1,2) may also be interpreted as 
dislocation densities. 

In the inclusion problems the input functions and the unknown 
functions are generally defined by [9] 

- - IT u mk > < k ■ '- 2) • (4) 

<i>l ( x ) - a 12 + - a ]2 , " °22 ' a 22 ’ ^ 

where u j and u m2 are the components of the displacement vector 
in the inclusion-free matrix at the location of the inclusion 
obtained from the given system of external loads, and a 12 , a 22 
are, respectively, the shear and the normal stress acting on the 
surface of the inclusion, the superscripts + and - again referring 
to the opposite faces of the inclusion. 

In (2) the kernels k^, (i,j = 1,2) are known functions 
which are bounded in plane problems and have logarithmic singu- 
larities in axisymmetric problems. In either case they are 
treated as simple Fredholm kernels. If the crack or the inclusion 
is completely imbedded in a homogeneous medium, then the constants 
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a.. (i,j = 1,2) vanish and the kernels in the third integral in 
J 

(2) are bounded. In this case the system of integral equations 

- 1/2 

has simple Cauchy kernels and the solution has r type singu- 

larity [2, 3, 5, 7]. If the imperfection lies along the interface, 
then the constants c.- k and d ijk are zero, the system of singular 
integral equations are of the second kind and the solution has 
the following typical oscillating singularity with a power -1/2: 

(^k ,a i j ) ~ r_1/2 sin [wlog(r/2a)] * (6) 

If an imperfection end terminates at an interface, from (2) 
it is seen that the kernels in the third integral become unbounded 
(as x' 1 ) as x and t approach the end point simultaneously. In 
this case the system of integral equations is said to have gener- 
alized Cauchy kernels, the power of the singularity of the solu- 
tion is no longer -1/2, and the solution of the system is not 
straightforward and requires much more care [73- 

In all the problems discussed here, it is not difficult to 
show that the index of the system of singular integral equations 
(2) is +1, hence the solution of the problem will contain two 
arbitrary constants [10]. These constants are determined by 
using the conditions 

a 

4> k (x)dx = 0 , (k = 1 ,2) . (7) 

-a 

In crack problems (7) expresses the conditions of si ngl e-val ued- 
ness of displacements, and in inclusion problems it constitutes 
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the conditions of static equilibrium of the inclusion. 

For the system of singular integral equations with simple 
Cauchy kernels the direct regularization technique in the type 
of problems under consideration becomes extremely cumbersome, 
and for the equations with generalized Cauchy kernels none exists. 
In practice these equations may be solved numerically in a very 
straightforward manner by using the methods described in [11-13]. 
In these techniques the function-theoretic method as described 
in [10] is applied directly to the integral equations and the 
fundamental functions of the system are obtained. Noting that 
these functions are the weights of certain Jacobi polynomials, 
and using the properties of the Jacobi polynomials, a set of 
Gauss-Jacobi type integration formulas are developed for the 
singular integrals and the system is solved numerically. Since 
the fundamental functions contain all the information regarding 
the singular behavior of the solution, and since they are deter- 
mined exactly, the resulting numerical solution preserves the 
correct nature of the singularity. 

3. CRACK PROBLEMS IN LAYERED COMPOSITES 

In plane or axisymmetric problems for layered composites 
containing a crack completely imbedded into a homogeneous medium 
and lying parallel to the interfaces, the fundamental functions 
and the solution of the system of singular integral equations 
are of the following form 
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5 


(8) 


w-|(x) = 


w 2 (x) 


= w ( x ) = (a‘ 


*V ,/2 


(J) k (x) = w(x) F k (x) , ( k = 1,2). 


(9) 


where the bounded functions F k (x), (k = 1,2) are determined from 
the solution of the integral equations. In this problem the 
increased number of layers simply affects the computation of the 
Fredholm kernels k,.(x,t), otherwise presents no conceptual dif- 

■ J 

ficulty. The stress intensity factors defined by 


k, (a ) = 1 i m [ 2 ( x - a)] 1/2 a (x,0) , 

x+a 


k ? ( a ) = 1 1m[2(x - a)] 1/2 o (x,0) , (10. a, b) 

L x-*a 

are easily shown to be related to the asymptotic values of -j and 
4> 2 as follows: 

k,(a) = " i 2 + 7 1 i m [ 2 ( a - x)] 1/2 ^ 2 ( x ) , 

• I ' ^ Y -V ^ 


k 9 ( a ) = - • 1 i m [ 2 ( a - x)] 1/2 ^(x) , (11. a, b) 

where y is the shear modulus, k = 3 - 4v for plane strain, and 
k = (3 - v)/(l + v) for generalized plane stress, v being the 
Poisson's ratio of the medium surrounding the crack. 

In this problem the strain energy release rate for the crack 
propagating in its own plane may be expressed as 

3U = ttQ ± _k! (k 2 + k 2 } . (12) 

3a 4y 1 2 

From the viewpoint of fracture propagation in brittle materials. 
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another important quantity is the "cleavage stress" o QQ given 


by [14] 

a 0e (r, e ) 


cosj [k, cos 2 f - f k 2 sine] + 0(r 1/r2 ) (13) 

(2 r)l/2 2 


where ( r , 0 ) are the polar coordinates at the crack tip. 

If the crack is on the interface, the system of singular 
integral equations is of the second kind and the fundamental 
function is given by [3] 


w ( x ) = (a - x) a (a + x) 3 , 
a = - j - iw , 3 = - j + iw , 


“ * H ’°9 ( H 


X) 

•\/ • 


Y 


(14) 


( y _ + K_y + ) - (y+ + K +0 
Y ” (y_ + K_y + ) + (y + + K + y_) 

where the subscripts + and - in the elastic constants refer to 
the materials on the y > 0 and y < 0 side of the crack, respec- 
tively. In this case, defining the stress intensity factors 


k-j and k 2 by 

k-, + ik ? = 1 i m [ 2 ( x 
1 c x+a 


a)] 


1/2 


(- 
' v 


x + a 


1W 

) [o yy (x.O) 


+ ia X y(x,0)] , 
(15) 


the strain energy release rate for a crack propagating along the 
interface may be expressed as [1] 


9U = TT 
3a 2 


1 + K 

~cy 


( k-j 2 + k 2 2 ) 




y+ - 

— i + ^ 

y_ + K_y + 


K + y _ - <_y + 
y + + K + y_" 


(16. a, b) 
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The functions f..(r,6) shown in (1) which give the angular 
variation of the stresses around the crack tip may be found, for 
example, in [14] for the homogeneous materials and in [15] for 
an interface crack. 

Figures 2-9 show some numerical results for various (plane 
strain) configurations involving the materials steel (E = 3x10 
psi, v = 0.3), aluminum (E = 10 7 psi , v =0.3) and epoxy (E = 
4.5xl0 5 psi, v = 0.35). The results are obtained for the per- 
turbation problem in which the crack surface tractions are 
assumed to be 

OyyCx.O) = - CT 0 . a xy (x ’ 0) = 0 * {17) 

Figure 2 shows that for the two (similar) half planes bonded 
through an epoxy layer which contains a symmetrically located 
crack, the stress intensity factors decrease with decreasing 
layer thickness. The opposite result is observed for the free 

layer, k-j and k ^ going to infinity as h 0 . The limiting value 

of the stress intensity factors, ii-|(l - v 2 ) / ( 1 " ^ or ^ ® 

is obtained by observing that the limit of the crack opening dis- 
placements in the layer 3 as h 0 is equal to that in the homo- 
geneous plane 2. If the crack is not in the mid-plane of the 

layer, k ^ is not zero. For a fixed ratio of (h/2a) - 1, the 
results are shown in Figure 3. In this and in the subsequent 
figures 0 is the probable cleavage angle obtained from [16] 

|j [(2r) 1/2 o 06 (r,e)] - 0 , 08) 

where the cleavage stress a Q0 is given by (13). The figure 
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clearly shows that, for this symmetric geometry, the propagating 
crack would tend to stay in the mid-plane of the layer. The 
figure also shows a measure of the strain energy release rate 
Wg defined by 


w 3 ■ — 2777-7 ( f?>, 

0 aircr ^(1 + Ko) 3 


k-| 2 + k 2 2 


(19) 


ac 


As h-j/h goes to 0 or 1 , the crack becomes an interface crack and 
the nature of the singularity changes. Thus at h-| = 0 and at 
h 1 = h, k-j and k 2 (considered as functions of h-j) are not con- 
tinuous. The quantity which is continuous is the strain energy 
release rate. Similar to (19) defining the measure of the, strain 
energy release rate for an interface crack (from (16)) by 


( k-j 2 + k 2 2 ) 


23 _ 


'23 


a o 2a 


2 c 2 3 U 3 / 9 U \ 

2/i . i ^9a%o 

ancr (1 + kJ 23 


( 20 ) 


the limiting values of W 3 may easily be found as follows: 


1 im W. 


W, 


h,*h ’ 3 " ’ 23 


1 i m W. 


W 


h r o ‘ 3 ■ 13 


( 21 . a ,b) 


where the constants c 23 and c ]3 are given by (16. b) and W 23 and 
W-. 3 are 1 
lem [4]. 


W are obtained from the solution of the interface crack prob- 

I o 


Results similar to that shown in Figures 2 and 3 are given 
in Figures 4 and 5 for three different materials. In this case 
too it is seen that (if u 3 < P-j > y 3 < for 3 given crack 

length as the layer thickness h decreases the resistance of the 
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bonded joint to fracture propagation increases (i.e., there is 
a reduction in W 3 and in the cleavage stress). Figure 6 shows 
the results for two bonded layers. Here, as the crack approaches 
the free surface, k-j , k 2 and become unbounded. For y 2 < » 

they all decrease as the crack approaches the interface, at 
h 3 = h 2 k-j and k 2 are again discontinuous, and W 2 = 2W 12 /c 12 
where W 12 may be evaluated from the solution of the interface 
crack problem (see Figure 9). 

Figures 7-9 show some of the results for the interface crack 
problem. Here k^- (°°) , (i = 1 ,2) and W 23 (°°) refer to the corre- 
sponding values for two bonded half planes with an interface 
crack (i.e., h = ») . In the case of three materials (Figure 7) 
if y 3 < u 2 , y 3 < y 1 , for a given flaw size again the fracture 
resistance of the bond increases with decreasing layer thickness. 
As seen from Figures 8 and 9, the opposite observation is valid 
if there is a free surface. Further results on the plane strain 
problem for a crack in layered materials may be found in [3,4]. 
The equivalent results for a penny-shaped crack are given in 
[5,6]. Similar results for an elastic fiber imbedded into an 
elastic matrix having perfect adhesion except for an axisymmetric 
interface crack are given in [8]. 

In the solutions given by the Figures 2-9 the adhesive layer 
is treated as an elastic continuum, and it is assumed that the 
flaw is a crack. In some problems this may be a realistic model. 
In some other problems, however, these assumptions would not only 
complicate the analysis, they may not even be realistic. This 
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would be the case, for example, in bonded materials in which the 
relative thickness of the adhesive is so small that the thickness 
variation of the stresses, in the adhesive layer may be neglected. 
In such problems the adhesive may be approximated by a combina- 
tion of shear and tension springs, which simplify the interface 
conditions considerably. Referring to the insert in Figure 10, 
for two dimensional problems these conditions may be expressed as 


- u l>/ h 3 = °xy /v 3 • 

(y = 

0) 


- v,)/h 3 = Oyy/E 3 , 

(y = 

o) 

(22 .a ,b) 


where h 3 , p 3 , E 3 are, respectively, the thickness, the shear 
modulus, and the Young's modulus of the adhesive, u^, v^ , (i 
1,2) are the displacements and 0 xy , c r yy are the tractions on the 
surfaces of the two materials bonded through the adhesive. 


For example, for the two half planes shown in Figure 10, 
defining the adhesive stresses by 

dyy(x,0) - p,(x) , o xy (x,0) = P 2 (x) (23 ' a>b) 

and using the solution for the elastic half plane [10], (22) may 
be expressed as 


/ P 2 (t) 1 °9 


dt + y / p -j ( t ) dt - bp 2 (x) + e - 0 


I P^t) log 


t - x 


dt - y / p 2 (t)dt - c p 1 (x) + 


/ p, (t) dt = P , I P 2 (t)dt = Q , 

-a “ a 


(24.a-d) 



where 


Y = Y-|/T2 ’ b ” h 3^ y 3 Y 2^ 5 C ” h 3^ E 3 Y 2^ ’ 

d = v 2 (0)/y 2 , e = u 2 (0)/y 2 > 


, K o - 1 K n "I 

t, ■ 1 < -Sr - ’ • Y: 


, K,+l K 9 + l 

]_ ( J + -A — ) , 

4tt v u 1 y 2 ' 


( 25 ) 


Pi ( X )> (i = 1,2) are the unknown functions and the unknown con- 
stants u 2 ( 0 ) and v 2 (0) are determined from the equilibrium con- 
ditions ( 24 . c ,d ) . The integral equations (24. a, b) have only 
logarithmic singularities; hence, the Fredholm theory is appli- 
cable and the functions p^, (i = 1,2) are bounded . 

Figure 10 shows the contact stresses for two elastic half 
planes bonded through an adhesive layer of width 2a and thickness 
hg* The figure also shows the results for the limiting case 
h = 0 for which the closed form solution is given by [15] 


Pl (x) 

p_ 

a + 1 1 

cos [wl og (|- 

+ i 
X X 
1 1 

2 tt 

„'/2 ( a 2 - x2) l/2 

XJ 
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X 

P_ 

a + 1 1 

si n [wl og (f- 

a 

+ 1 

X X 

1 1 

V* 

2 TV 

a 1 ' 2 (a 2 - x 2 ) 1/2 


where 

a = ( <1 u 2 + y-j ) / ( < 2 V1 1 + y 2^ * w = h\ 109 “ ‘ 

It should be noted that, for the configuration given in 

*From (24) it can be shown that at x = + a even though Pl -(x), 

(i = 1,2) is bounded, the derivative dp^dx has a logarithmic 
singularity, i.e., d Pi /dx « K log(a-x) for small values of (a-x). 



Figure 10 if the relative thickness of the adhesive, h 3 /a is 
increased, the contact stresses become more "uniform" and their 
peak values decrease. Table 1 shows the values of the stress 
concentration factors" for three different adhesive thicknesses. 
Further application of this concept to the stress analysis of 
the adhesive joints may be found in [17]. 

Table 1. Stress concentration factors 
for two bonded half planes. 

h/a 0.004 0.010 0.040 

a yy (a,0)/(P/2a) 4.240 2.900 1.672 

a (a,0)/(P/2a) -0.7020 -0.3322 -0.1148 

xy 

4. A CRACK TERMINATING AT AND GOING THROUGH AN INTERFACE 

The solution of the problem of a crack propagating toward 
and terminating at a bi material interface was given in [7]. 

From the viewpoint of fracture propagation, in this problem the 
quantities of interest are the stress intensity factors, the 
power of the singularity, a (see equation (1)), and the angular 
distribution of the stresses around the crack tips. For the 
crack tip imbedded in the homogeneous medium a = 1/2 and the 
6-dependence of stresses are well-known (see, for example, [14]). 
In [7] a was obtained directly from the integral equations by 
using the f uncti on-theoreti c approach. It may also be obtained 
by using the technique of the eigenfunction expansion [18]. 

The 6-dependence of the stresses around the crack tip ter- 

minating at the interface is dependent on the elastic constants 

u K (i = 12 ) and is obtained from the solution of the 
M i ’ i 

-17- 



problem [7]. Figures 11-13 show the results for two bonded semi 
infinite planes where material 1 contains a semi - i nf i ni te crack 
perpendicular to and terminating at the interface. The plane is 
loaded by symmetric wedge forces P applied to the crack surfaces 
at a distance r Q from the crack tip. In this problem for small 
values of ( r / r 0 ) the stresses may be expressed as 


a n (r,e) ■ E- (^) G,j(e) + 0 [(r/r o ) 2 ] . 

(0 < a < 1, (1 ,j) = (r,0), 0 < 9. < ", Re l>2 ] > 0). (27) 
The functions G,, are shown in Figures 11-13 where m = V 2 /' y l and 

• J 

the corresponding material combinations are given in Table 2. 

The table also shows the stress intensity factor k defined by 

k = lim n r a a 206 (r,O) (28) 

r->0 

which can be evaluated in closed form [7] and which, in the homo 

i 

geneous medium, is' given by 


k = P (— ) 

* ^o 


1/2 


(29) 


Table 2. Stress intensity factors 
for semi’-i nf i ni te crack. 


Material m = 1*2^1 k/ ( P/vr Q ) 


Boron-Epoxy 
A1 umi num-Epoxy 
A1 umi num-Al uminum 
Epoxy-Al uminum 
Epoxy-Boron 
Epoxy-Ri gi d 


0.0072 

0.0511 

0.0433 

0.1197 

1 .00 

0.4502 

23.077 

0.8403 

1 38.462 

0.8905 

oo 

0.9015 
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In the general case (for symmetric loading) the asymptotic values 
of the stresses for small r may be expressed as 


a,,(r,e) = f 1 j ( 0 ) + 0(r“ 2 ) , 

ij /2 r J 

( Rea 2 > 0 , i , j = r ,6) , (30) 

where the functions f.j are related to G^. shown in Figures 11-13 
through 

f . . (0 ) = G..(6)/G (0) , { i » J = r,0). (31) 

I vj * J ' 


In the case of a finite crack perpendicular to the interface, 
the stress intensity factors are shown in Figure 14. Here the 
external load is a uniform pressure, p Q applied to the crack 
surface and the two materials considered are Aluminum (v = 0.3) 
and Epoxy (v = 0.35) (see Table 2 for modulus ratio). The power 
of the singularity, a for the crack tip terminating at the inter- 
face (i.e., c = a ) is given in Tabled 3. Note that, because of 
the variation in k, there is a slight difference between the 


values of a corresponding to plane strain and plane stress. 


Table 3. Stress intensity factors and the power of singularity 
for a crack terminating at the interface. 


V 1 *! 


Plane Strain 


Plane Stress 


k(a) ! k (b) 
Pc/^o 


k(a ) k(b) 

Po^ Po^ 


23.077 


0.3381 

0.8258 


2.6237 

0.0700 


0.8827 
1 .3552 


0.2890 4.1760 0.8789 

0.8230 0.0744 1 .3525 


0.0433 
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For this case Table 3 also shows the stress intensity factors at 
the crack tips r = a = 0 and r = b which are defined by 

k(a) = lim /2 r a a 200 (r,O) , 
r+0 

k(b) = lim /2 ( r - b) o-| Q0 (r , it ) . (32. a, b) 

r+b 


Due to the deviation in a from 1/2 as given in Table 3, the limits 
of the curves k ( a ) shown in Figure 14 as c -> a Q are 0 and °° for 
y 1 < y 2 and y-j > u 2 , respectively. 


To solve the problem of a crack going through the interface, 
first the problem of two cracks (one in each material) perpendic- 
ular to the interface is considered. The problem is again for- 
mulated by using the Mellin transform. In the resulting system 
of singular integral equations the input functions are the crack 
surface tractions. For the case in which the two cracks are 
imbedded into the adjacent homogeneous media, the kernels of the 
integral equations have only Cauchy type singularities. Hence 
the solution of the problem is rather straightforward. When the 
inner crack tips go to the interface, the cracks join and become 
a through crack. In this limiting case, even though the integral 
equations remain unchanged, the singularities of the kernels 
become of the generalized Cauchy type. This means that r - 0 

(as well as r = b, and r = bj is an irregular point around which 

* 

the solution will have a singularity of the form r a , (0<a<l/2) . 


^ T h i s is the case for most of the practical material combinations 
However, for some unusual combinations the stresses at r = 0 may 
be bounded or may have only a logarithmic singularity (see L 1 9 J 
for a thorough discussion). 
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The system of characteristic equations giving a and the powers 
of the singularity at the two (outside) ends of the crack is 
obtained by applying the f uncti on- theoreti c method directly to 
the system of singular integral equations. 


In the examples considered in this paper, it is assumed 
that the material is subjected to tensile loads parallel to the 
interface and sufficiently away from the location of the crack. 
For the perturbation problem it may then be assumed that the 
crack surface tractions are uniform pressures related by 

p,/p 2 - e,/e 2 (33) 

TT 3 7T 

where the subscripts 1 and 2 refer to materials 1 (^ < 9 < 2 ~) 
and 2 (- | < 9 < |) » respectively. 

The stress intensity factors for two cracks imbedded into 
the adjacent bonded planes are shown in Figure 15. In this 
example Material 1 is aluminum, Material 2 is Epoxy, the crack 
lengths 2^ = b, -a,, 2* 2 = b 2 - a 2 are equal and fixed, a 2 is 

constant, a 2 = 2* 2 is the length unit, and the variable is the 
distance of the first crack to the interface, or (b 1 + z ] )/2* 2 - 
In the special case of the homogeneous medium containing two 
col 1 inear cracks of equal lengths and subjected to a uniform 
pressure p Q on the crack surfaces, the stress intensity factors 

are given by [20] 


k(a) = 
Po^ 


— ^ TP? 5 

(b - a ) [a ( b + a ) ] ' ^ 
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(34. a ,b) 


k(b) /? b 2 ri - E ( m ) / K ( m ) ] 

p Q /I (b - a ) [ b( b + a ) ] ' ^ 

(2£ = b - a , m = 1 - a 2 /b 2 ) , 

where 2a and 2b are the distances between inner and outer tips 
of the cracks, respectively. For example, for b = 2a we obtain 

k( a ) = 1.01762 p Q £ 1/2 , k ( b ) = 1.01249 P Q £ 1/2 • 


Figures 16-22 show some of the results obtained for the 
material combination Aluminum-Epoxy. The stress intensity fac 
tors shown in these figures are defined as follows: 


k, = lim /2(r - b ] ) a 100 (r,TT) , 

r+b-j 

k 2 = lim /2 ( r - ’ 

L r->b 2 

k = lim r“ o lee (r,j) . 
r+0 


k = lim r 
r r+0 


a 


a l r6^ r ’2^ 


( 35 . a-d ) 


The factors k 0 and k r are a measure of the contact stresses 
around r » 0 and may be used in studying the question of debond- 
ing crack initiation. The stress intensity factors shown in the 
figures are normalized with respect to p-| & 
where 2SL is the total crack length. 


, 1/2 


p 2 &^ 2 , and p-|& a , 


In Figures 16-18 material 1 is aluminum, material 2 is epoxy, 
the penetration of the crack into 1 is fixed and is taken as the 
length unit, i.e., b! = 1, and the crack length 2£ (or the crack 
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penetration into the second medium) is the variable. For plane 
strain case and for the material pair under consideration the 
power of singularity at r = 0 (i.e., for two bonded quarter 
planes) is found to be a = 0.273692. In the limiting case of 
l = 0.5, b 2 = 0 and the problem reduces to that of a crack ter- 
minating at the interface for which the results are given in 
Table 3. Since at r = 0 the power of singularity for this limit- 
ing case (i.e., a = 0.8258) is greater than 1/2, as seen from 
Figures 16 and 17, the stress intensity factors k 2 , k Q , and k r 
become unbounded as £ - 0.5. For the same material combination 
the crack surface displacement is shown in Figure 18. Here the 
functions u l0 (r,ir), (0 < r < b } ) and u 2Q (r,0), (0 < r < b 2 ) are 
plotted against the dimensionless variables and defined by 

C , . 2^-1 . < 2 -^-’ • (36 - a - b) 

As seen from the figure, because of the singularity at r = 0, it 
turns out that the displacement has an unbounded derivative at 
r = 0 as well as at r = b^ , (i = 1,2). 

Figures 19 and 20 show the results analogous to that given 
by Figures 16 and 17, the difference being that here material 1 
is epoxy and material 2 is aluminum. In the limiting case of 
l = 1/2 the power of the singularity is a = 0.3381 (see Table 3) 
which is less than 1/2 and is greater than the power of the 
singularity a = 0.273692 at r = 0 for £ •> 1/2. Thus, as seen 
from Figures 19 and 20, as £ goes to 1/2, k 2 goes to zero and 

k and k Q become unbounded. 

r t) 
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Further results for the material Combination Aluminum (1) - 
Epoxy (2) are shown in Figures 21 and 22. In this case the total 
crack length 2£ is constant and the crack eccentricity c is 
selected as the variable. 

5 . PENNY-SHAPED CRACK IN A FILAMENT-REINFORCED ELASTIC MAT RI_X 

In filament-reinforced composites the spacing and the orien- 
tation of the filaments are generally random. However, in some 
cases it is desirable and, to a certain extent, possible to con- 
trol the orientation of the filaments in the matrix. In this 
Section we will describe some of the results of a study on the 
stress distribution in a filament-reinforced elastic matrix con- 
taining a penny-shaped crack [21]. It will be assumed that 
(a) the matrix is "sparsely" reinforced by a finite number of 
filaments which are symmetrically distributed around the crack 
and are oriented perpendicular to the plane of the crack, (b) 
z = 0 is a plane of symmetry for the external loads and for the 
geometry, and (c) the filament radius r Q is small compared to 
its length 2c. Under these assumptions the problem can be for- 
mulated in terms of an integral equation with essentially the 
interface shear stress between the filament and the matrix as 
the unknown function. The main interest in this problem is in 
the evaluation of (a) the stress intensity factor at the crack 
periphery as a function of 6, (0 < 0 < 2tt), (b) the contact 
stress between the filament and the matrix as a function of z, 
(-c < z < c), and (c) the stress concentration factor in the 

f i 1 ament . 
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Figures 23-25 show the 0-dependence of the stress intensity 

factor, K-j . In the examples given in this Section it is assumed 

that the medium is subjected to uniaxial tension in z-direction, 

a - o away from the location of the crack and the filaments, 
zz o’ J 

Figure 23 shows the variation of the stress intensity factor in 
the case of a single filament for various values of filament-to- 
matrix modulus ratio, E^r/ E and for f i 1 ament-to-crack center dis- 
tance, b. It is seen that the effect of the filament on the 
stress intensity factor becomes more significant as the distance 
b decreases and as the modulus ratio E^/E increases. 

The effect of the number L of symmetrically located (iden- 
tical) filaments on the stress intensity factor for a given modu- 
lus ratio, ( E f /E) = 100 is shown in Figure 24. The figure indi- 
cates that for b i 1.5 and L z 6 the 6-dependence of the stress 
intensity factor may be neglected. Figure 25 shows the effect 
of the modulus ratio, E^/E on K-| for L = 2 and L = 4. 

Finally, the effect of the modulus ratio on the maximum 

tensile stress in the filament is shown in Figure 26. Here 

~o f (0) is the stress in the filament at z = 0 and a Q is the 
c 

stress in the matrix away from the crack-filament region. For 
large values of the filament length 2c and the distance b we 
have (a c f (0)/c o ) * ( E f /E) . For comparison the figure also shows 
the stress at the same location in the absence of any filaments 
(i.e., for E f = E) obtained from the axisymmetric penny-shaped 
crack solution [22]. Further results on this problem may be 
found in [21]. An approximate solution of a special case of the 
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f / 

problem is given in [23] where the 6-dependence of the stress 
state in the composite medium is neglected. 


6. INCLUSION PROBLEMS IN BONDED DISSIMILAR MATERIALS 

In some composites material imperfection may be in the form 
of an elastic inclusion rather than a crack. In this case it 
may be expected that due to the high stress concentration, the 
neighborhood of inclusion borders would be the likely location 
at which the fracture nucleation would occur. To give some idea 
about the nature of the problem, in this Section we will discuss 
the solution of an idealized plane inclusion problem. It will 
be assumed that the inclusion is oriented parallel to the inter- 
face of two bonded elastic half planes and the thickness of the 
inclusion is very small compared to its lateral dimensions. 

Thus, analytically, it can be approximated by a singular surface 
across which the displacement vector is continuous and the stress 
vector suffers a discontinuity. The formulation of the problem 
then leads to a system of singular integral equations. 


In the case of a perfectly rigid inclusion, the solution of 
the plane problem gives the stress state \in the close neighbor- 
hood of an end point as follows (see [9]): 


°yy (r,9) 


[tti (k i co 4 + r-H k 2 s1 4' + 0(rl/2) ’ 


(2r) 


o ( r , 0 ) 
xy v 


a ( r , 0 ) 
xx v 


V? (- f-TT k l si 4 + k 2 C0 4> + °< rl/2 > • 


(2r) 


T72 


7172 <- H-T k l c °4 + Kt k 2 si 4> + 0(rl/2) 


( 2 r ) 


1/2 
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re i0 = x - a + iy , (-ir < 6 < tt) , (37.a-c) 

where a is the half-length of the inclusion. If the distance, 
h from the inclusion to the interface becomes zero (i.e., if the 
inclusion is located on the interface), the problem can be solved 
in closed form [9]. In this case, similar to the interface crack 
problems, the stresses at the end point have the well-known 
oscillating singularity; hence, at h = 0, the stress intensity 
factors k-j and k 2 (considered as functions of h) are discontinuous. 

Perhaps a more realistic model for the inclusion would be an 
elastic layer with known extensional stiffness. Since the rela- 
tive thickness of the inclusion is assumed to be very small, its 
bending stiffness will also be very small and may be neglected. 

For this "elastic membrane" type of inclusion the stress state 

i 

around the end point may be expressed as [9] 


o ( r ,0 ) = 

yy 


° xy (r,e) 


a (r ,0) = 
xx v 


( 2 r ) 


1/2 


( 2 r ) 


T72 


cos|- + 0(r^ 2 ) , 


sin| + 0(r 1/2 ) , 


( 2 r ) 


3 k -j -j ~2 cos| + 0(r 1/2 ) 


(38.a-c) 


Here as h •> 0, unlike the problem of rigid inclusion, the inte- 
gral equation does not change its character and remains as a 
singular integral equation of the first kind. Hence the singular 
behavior of the stress state maintains its simple r character 

( see Fi gure 28) . 
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Figures 27-29 show some of the results obtained for rigid 
and elastic inclusion. In these examples it is assumed that the 
medium is subjected to symmetric loading away from the location 
of the inclusion. Thus the displacements in the medium without 
the inclusion and at the same location as the inclusion may be 
expressed as 

«„ * - v . % ' ° • (39) 

The constant k Q shown in the figures which is used to normalize 
the stress intensity factors is defined by 

k 0 - H 2 V 1/2 (k 2 ■ 1)/k 2 ’ <40) 

where u 2 and * 2 are the elastic constants of the medium surround- 
ing the inclusion. 

Figure 27 shows the results for a rigid inclusion in two 
bonded elastic half planes and that in a half plane. Similar 
results for an elastic inclusion are shown in Figure 28. Figure 
28 also shows the stress intensity factor for an inextensible 
inclusion with zero bending stiffness. The values of k for 
h + o are obtained from the interface inclusion problem or from 
the problem of a stiffener bonded to an elastic half plane [24]. 
The figure indicates that the stress intensity factor is highly 
dependent on the stiffness of the inclusion. 

In this problem too one could determine the orientation of 
the probable fracture nucleation by considering the cleavage 
stress o aR as a function of 9. For example, for the elastic 
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inclusion problem from (38) we obtain 


a., (r.e) - K-T 2 (2cos§ COS20 - cosf - sinf sin26) 

ee ( 2r ) ' c 

+ olr 1 ^ 2 ) , (-ir < 0 < ir) . (41 ) 

From Oo qq / 30) =0 we find 6 = 0°, 74.3°, 151°. The angular 
variation of a Q0 is shown in Figure 29. It is seen that if the 
medium is loaded perpendicular to the inclusion the fracture 
initiation will most likely be collinear with the inclusion with 
possible secondary radial cracks along 6 = ± 151°. On the other 
hand if the medium is loaded parallel to the inclusion the frac- 
ture would most likely nucleate along (roughly) 0 = ± 74.3° 
radial lines. 
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Crack in the layer bonding two different half planes 
(h = 2a, Materials: same as in Fig. 4.) 





1.2 
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Fig. 7. Interface crack in two half planes bonded through an 
elastic layer, k-j , V >2 • stress intensity factors, W 22 
strain energy release rate. (Materials: 1 : Steel , 

2: Aluminum, 3: Epoxy; load: crack surface pressure.) 









h 2 /2a 


Fig. 9. Interface crack between two bonded elastic layers. 

(h = 3h 2 ; Materials: 2: Steel, 3: Aluminum; load: 
crack surface pressure.) 
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Fig. 10. 


Contact stresses for two elastic half planes bonded 
through a thin adhesive layer. 





Fig. 12. Angular variation of o r g around the crack tip termi 
nating at an interface. 
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Fig. 13. Angular variation of a rr around the crack tip termi 
nating at an interface. . 
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Fig. 14. 


Stress intensity factors for a crack perpendicular to 
an interface. (Materials: Aluminum and Epoxy; load: 


crack surface pressure. 


) 
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Fig. 15. Stress intensity factors for two collinear cracks 

imbedded in the bonded adjacent materials. (Materials 
1: Aluminum, 2: Epoxy; a 2 = 1 = 2i 2 = 251-, , b 2 = 2a 2 ; 
load: crack surface pressures p-| and p 2 : ( p -, / P 2 ) = 
(E,/E ? ).) 









Fig. 17. Stress intensity factors for cleavage and shear stress 
on the interface for a crack going through the inter- 
face. (Materials: 1: Aluminum, 2: Epoxy, b ] = 1, 

% = ( bi + b 2 ) / 2 ; load: crack surface pressures P 1 , P 2 - 





Crack surface displacement for a crack going through 
the interface. (Materials: 1 : Aluminum, 2. Epoxy, 







Stress intensity factors for cleavage and shear stresses 
on the interface. (Materials: 1: Epoxy, 2: Aluminum, 
l = (b, + b~)/2, load: crack surface pressures p 1 , P 2 *) 



Stress intensity factors for a crack going through 
the interface. (Materials: 1: Aluminum, 2: Epoxy, 
21 = b-j + b 2 = constant, c = | b -j - b 2 |; load: crac 
surface pressures.) 








Fig. 23. Stress intensity factor on the border of a penny-shaped 
crack in a matrix reinforced by a single filament. 
(E.p/E: f i 1 ament-to-matrix modulus ratio, b: filament 
distance to the crack center, crack radius = unity, 

2c: filament length, r Q : filament radius, 0: angle 
measured from r axis intersecting the filament; load: 
uniaxial stress o n at infinity.) 




QL/tt 

Fig. 24. The effect of filament distance b, and the number of 

filaments L on the stress intensity factor for a penny- 
shaped crack in a filament-reinforced matrix. 
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Fig. 25. The effect of modulus ratio E f /E on the stress intensity 
factor for a penny-shaped crack in a filament-reinforced 
matrix. 
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Fig. 26. The stress in t 
tion of f i 1 amen 
shaped crack in 


E,‘/E=5 


Ef/E = l 
No Filament 

I L_ 

2.0 3.0 

b 

e mid-plane of the filament as a func- 
distance b. (Radius of the penny- 
the matrix is the length unit.) 





Fig. 27. Stress intensity factors for a flat rigid inclusion 
parallel to the bi-material interface. (Materials: 
1: Aluminum, 2: Epoxy; load: displacement in the 
inclusion-free medium 2, u Q = - e Q x, v Q =0; 

k o = y 2 e o al/2 ^ K 2 “ 1 ) /k 2*^ 5^ 
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Fig. 28. Stress intensity factors for a flat elastic and for a 
flat i nextens i bl e inclusion lying parallel to the 
interface. (Materials: 1: Aluminum, 2: Epoxy, load: 

= P 2 e o al/Z ^ ,c 2 " 1 K 2 * ^ 
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